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This paper establishes expectation and variance asymptotics for statistics of the Poisson-Voronoi 
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1. Main results 

The Poisson-Voronoi mosaic is a classical and prominent example of a random mosaic 
and is used in a wide range of fields, including astronomy, biology, material sciences and 
telecommunications. If V\ is a Poisson point process on Q := [—1/2,1/2] d whose intensity 
measure has density Ak(-) with respect to the Lebesgue measure (d > 2, A £ (0, oo) and 
n is a continuous function on Q bounded away from zero and infinity), the Voronoi cell 
v(x ) := v(x,P\) associated with x GV\ is the set of all z £ Q such that the distance 
between z and x is less than the distance between z and any other point of V \. Clearly, 
v{x) is a random convex polytope and the collection of all v{x) with x GV\ partitions Q 
and is called the Poisson-Voronoi mosaic of Q. 

Let A C Q be a full-dimensional admissible set whose boundary has positive and finite 
(d— l)-dimensional Hausdorff measure. Admissible sets, formally defined in Section 2, 
include in particular, convex sets, sets of positive reach, differentiable manifolds with 
smooth boundary as well as certain finite unions of such sets. Given such A C Q, the 
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Poisson-Voronoi approximation PVa(A) of A is the union of all Voronoi cells v(x) with 
x € A, that is, 


PVa(^) := (J v(x). 


xGVxnA 


Typically, A is an unknown set having unknown geometric characteristics such as volume 
or surface area. Notice that PVa(AI) is a random polyhedral approximation of A, which 
closely approximates A as A becomes large. One might expect that the volume and surface 
area of PVa(AI), respectively denoted by V\(A) and S\ (A), also closely approximate the 
volume and surface area of A. Our first goal is to show that this is indeed the case, 
though the surface area asymptotics involve a universal correction factor, denoted by C 2 
in the sequel. For sets A which are convex or which have a smooth boundary, first-order 
asymptotics have been previously established in [9, 21, 25, 31]; second-order asymptotics 
for sets A having a smooth boundary are given in [31], while [28] provides second-order 
inequalities when A is a convex set. We extend the limit theory of these papers and 
obtain first- and second-order asymptotics whenever A belongs to the more general class 
of admissible sets. In particular, we show that the variance asymptotics for Va(^) are 
proportional to the K-weighted surface content of A , resolving a conjecture implicit in 
Remark 2.2 of [25]. The approach relies on a general and far-reaching Steiner-type formula 
from [10], together with stabilization properties of geometric functionals of the Poisson- 
Voronoi mosaic. 

In the sequel, we write /(A) ~ cg( A) for real-valued functions / and g and constants 
c G [0, oo ) if liniA^oc /(A)/g(A) = c. Throughout, we denote the s-dimensional Hausdorff 
measure by 'H s , s G [0,oo). Furthermore, we say that dA contains a subset T of differ¬ 
entiability class C 2 with / H d ~ 1 (T) G (0, oo) if T C dA is an open and twice differentiable 
(d— 1 [-dimensional sub-manifold in R d in the usual sense of differential geometry. Finally, 
for 7 G R we define the n-weighted surface content 



Observe that H. d .}{dA) reduces to the usual surface content 'H d x (dA) of dA if either 
7 = d and n is arbitrary or k = 1 and 7 G R is arbitrary. 

Theorem 1.1. There are constants ci,C 2 G (0,oo) depending only on the dimension d 
such that 


EV x (A)-V{A)~ Cl \- 1 / d H d - 1 {dA) and ES X (A) ~ c 2 U d J_ 1 (dA). 


Moreover, there are constants 03 , 044 , 04 ^ G [0, 00 ) depending only on d such that 


Var[R A (^)] ~ c 3 \- 1 - 1/d U d - 1 (dA) 


and 
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If dA contains a subset T of differentiability class C 2 with 'H d 1 (r) £ (0,oo), and if 
k= 1, then C 3 and C 4 := c^i + c^i are strictly positive. 

Next, we turn to other metric parameters of the Poisson Voronoi approximation, which 
can be handled by our general set-up. To this end, for I £ {0,..., d — 1} denote by 
skeb(PVA(A)) the union of all ^-dimensional faces belonging to <9 (PVa(M)), the boundary 
of PVa(M), and let H^\A) be the ^-dimensional Hausdorff measure of skeb(PV,\(^l)). 
More formally, if Pe(P) stands for the collection of ^-dimensional faces of a polytope P, 
then 

£ «'</>■ 

xGVx feJ 7 e(v(x)) 

X ^ A / C 0(pVa(A)) 

Note that H^ 1 ^(A) coincides with S\(A) considered in Theorem 1.1. 

Theorem 1.2. Let I £ {0,..., d— 1}. Then there are constants C 5 £ (0, 00 ) and C 6 ,i, ce.i £ 
[ 0 , 00 ) depending only on d and t such that 

E H^\A) ~ c h \ x - xld - lld U d -l(dA) 

and 

Var[Pf (A)] ~ X^^-^ice.iH^dA) + c 6 , 2 P^(£>A)). 

If dA contains a subset T of differentiability class C 2 with 'H d ~ x {T) £ (0, 00 ), and if 
n=l, then cq := C 6 ,i + C 6, 2 is strictly positive. 

With the exception of the number of vertices on 5 (PVa(A)), we have investi¬ 

gated only metric parameters of the Poisson-Voronoi approximation, namely the volume, 
the surface area and the Hausdorff measure of lower-dimensional skeletons. On the other 
hand, the combinatorial complexity of PVa(A) is also of interest. For example, it is nat¬ 
ural to ask how many ^-dimensional faces (£ £ {0,..., d — 1}) belong to <9 (PVa(A)). In 
contrast to volume and surface area, combinatorial parameters of the Poisson-Voronoi 
approximation have apparently not been studied in the literature. The general theory 
developed in Section 2 allows us to investigate such parameters. To state the result, 
for £ £ {0,... ,d — 1} we let f^’ 1 (A) be the number of ^-dimensional faces belonging to 
<9 (PV a (A)). Note that /f } (A) = P{ 0 ) (A). 

Theorem 1.3. Let £ £ {0,..., d— 1}. Then there are constants C 7 £ (0, 00 ) and esp, cs , 2 £ 
[ 0 , 00 ) depending only on the dimension d and on £ such that 

Efi e \A)^c 7 X 1 - 1 / d n d - 0 1 (dA) 

and 

Var[./f {A)] ~ X'-WfaUifoHdA) + c 8 , 2 H d ^ 0 (dA)). 
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If dA contains a subset T of differentiability class C 2 with Hr 1 (r) £ (0,oo), and if 
k = 1, then c$ := Cs,i + cs^ is strictly positive. 

Next, we consider certain functionals of Voronoi cells intersecting only a part of the 
boundary of A. Formally, given an admissible set A and Ag C dA such that H d ~ 1 (Ag) £ 
(0, oo), define the Poissson-Voronoi zone PVZ a (A 0 ) of Ag by 

PVZ a (A 0 ) := U v{x). 

xeV\ 

u(x)nAo^0 


w p\ 

Given £g {0, 1}, let f^'(Ao) denote the number of ^-dimensional faces of 

PVZa(Ao). We emphasize that this construction is very similar to the construction of 
a zone in a hyperplane arrangement; see [16]. Following these classical ideas, we define 
the complexity of PVZ a (A 0 ) as Co a (A 0 ) := f^\Ag) + ••• + /} d_ 1 ^(Ao). The zone the¬ 
orem in discrete geometry (see Theorem 6.4.1 in [16]) asserts that the complexity of a 
zone of an arbitrary hyperplane arrangement is of surface-order. Our next result shows a 
similar surface-order behaviour for the expectation and the variance in case of a random 
Poisson-Voronoi zone. 

Theorem 1.4. There are constants eg £ (0,oo) and cio.i,cio ,2 £ [0, oo) depending only 
on d such that 

EGo x (A 0 )^cg\ 1 -^ d H d - 0 1 (A 0 ) 

and 

Var[Co A (A 0 )] ~ A 1 ^(ciop'H^^Ao) + cig^'H d K ffj ) {A 0 )). 

If Ag contains a subset T of differentiability class C 2 with 'H d ~ 1 (Y) £ (0,oo), andifn = 1, 
then cio := cio,i + cio ,2 is strictly positive. 

Another application of our results concerns the iterated Poisson-Voronoi approxima¬ 
tion, defined recursively as follows: 

PVl 1 } (A):=PV A (A) and PV^(A) := PV^PV^-^A)) 

for integers n > 1 (note that the intensity used in the nth iteration is nA, where A > 0 is 
fixed). By V A "\ S' A n) and we denote the volume, the surface area and the number 

of ^-dimensional faces (£ £ {0,..., d— 1}) of the nth iterated Poisson-Voronoi approxima- 
tion, respectively. Moreover, by Hf we indicate the ^-dimensional Hausdorff measure 
of the ^-skeleton of PV^(A), £ £ {0,... ,d}. Note that our construction of the iterated 
Poisson-Voronoi approximation is close to that of so-called aggregate mosaics introduced 
in [29]. The expectation analysis of functionals of the iterated Poisson-Voronoi mosaic 
yields the following result. Variance asymptotics are less tractable and we shall omit 
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them. For simplicity, we shall assume that the Poisson point process V\ is homogeneous 
with k = 1 . 


Theorem 1.5. Suppose that k= 1 and let Ci and c 2 be the constants from Theorem 1.1, 
C 5 the constant from Theorem 1.2, and cj the constant from Theorem 1.3. Put c 2j „ ~ 
1 + C 2 + c| H + df ~ 1 for integers n > 1. Then 


EV A (n) 

- V(A) r 

' J cic 2 ,nA i/ d U d 1 (dA), 

E S<" : 

'~S(A) r 

- C2C2, n X~ 1/d n d ~d- 1 (dA), 


EH £ f {n) r- 

- c b C2, n \ l - lld - l/d U d -l{dA) 


Ef an) r 

- C 7C2,nX 1 - 1/d -H d - 0 1 (dA). 


Remarks. 

(i) Theorem 1.1 (related work). The set PV,\(A) was introduced in [13] where it was 
shown that limA-j-oo Vol(AAy4>,) =0 almost surely, but only when d= 1. This almost 
sure limit was extended in [21] to all dimensions d > 1. When V\ denotes a homogeneous 
Poisson point process on R d having intensity A, we have that V\(A) is an unbiased 
estimator of V(A) (cf. [25]), which makes PVa(AL) of interest in image analysis, non- 
parametric statistics and quantization; see also Section 1 of [13] and Section 1 of [9]. 

(ii) Invariance of limits with respect to geometry. The common thread linking our 
results is that the first- and second-order asymptotic behaviour of our functionals are 
geometry independent. By this we mean that the mean and variance asymptotics are 
not influenced by the precise geometric structure of the given admissible set A , but are 
rather controlled only by the k -weighted surface content of A. 

(iii) The constants in Theorems 1.1-1.5. The explicit dependency of the constants 
Ci,i > 1, in Theorems 1.1-1.5 on the dimension d and the parameter I is given explicitly 
in the general results of Section 2, especially the upcoming limits (2.16) and (2.17). More 
precisely, let 'P] lom be a homogeneous Poisson point process on R d of unit intensity and 
put R^. -1 :=R d_1 x R + . Let 


PV(Rt 1 ) := U «(*) 

^.g-phompjjd-l 

be the Poisson-Voronoi approximation of Rl . Then the general results show that the 
expectation and variance asymptotics are controlled by the K-weighted surface content 
of A as well as by the expected behaviour of metric and combinatorial parameters of the 
simpler object PV(R(Jr 1 ). Finding explicit numerical values for the constants Ci,i> 1, 
arising in expectation and variance asymptotics is a separate problem which we do not 
tackle here. 
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(iv) Extensions of Theorems 1.1-1.5. By Theorem 2.1 below, the expectation asymp¬ 
totics in Theorems 1.1—1.5 may be upgraded to a weak law of large numbers holding in 
the L 1 - and L 2 -sense. 

(v) General surface-order results. Although Theorems 1.1—1.5 only deal with statis¬ 
tics of the Poisson-Voronoi approximation, we emphasize that they follow from general 
theorems (presented in Section 2 below) for general surface-order stabilizing function¬ 
als. These general theorems are applicable in a wider context, establishing, for example, 
expectation and variance asymptotics for the number of maximal points in a random 
sample, as described in Remark (iii) after Theorem 2.2. 

The rest of this paper is structured as follows. In Section 2, we make precise our 
framework, in particular, we introduce the class of admissible sets and score functions. 
We also state there two general theorems which yield Theorems 1.1-1.5. Their proofs 
form the content of Section 3, while Section 4 contains the proofs of Theorems 1.1-1.5. 
Section 5 establishes the asserted variance lower bounds in Theorems 1.1-1.4. 


2. Framework and general theorems 

Let V\ denote a Poisson point process on R d for some d > 2 whose intensity measure 
has density A n with respect to the Lebesgue measure on R d , where A € (0,oo) but now n 
is a bounded function on R d not necessarily bounded away from zero. Furthermore, let 
A C l d be a closed set such that its boundary dA has finite (d— l)-dimensional Hausdorff 
measure. We consider in this section general statistics of the form 

fa,P x ,dA), (2.1) 


where £ is a certain score function, which associates to a point x G V\ a real number, 
which is allowed to depend on the surrounding point configuration V\ as well as on the 
set A via its boundary dA. To introduce a re-scaled version and to simplify notation, we 
use the abbreviation ,dA) \= ( t {\ 1 / d x,\ 1 ^ d V\ 1 X l ^ d {dA)) and define 

H*V> x ,dA):= Ux,Vx,dA). ( 2 . 2 ) 

x&Vx 

The focus of this paper is on score functions which depend on the geometry of the set 
A in that t;(x,V\,dA) decays with the distance of x to dA. Moreover, we require £ to 
satisfy a weak spatial dependency condition. 

To make the framework precise, we first introduce terminology, including the collection 
A (d) of admissible sets A C ! d as well as the collection H of admissible score functions. 
The reader may wonder about our choice of admissible sets. The admissible sets described 
below have the attractive feature that their so-called extended support measures are ‘well- 
behaved’ and satisfy a Steiner-type formula (2.3), which is a far reaching consequence of 
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the classic Steiner formula. This key formula, proved in [10], essentially replaces the co¬ 
area formula applicable in the surface-order asymptotics of functionals of sets A having 
a smooth boundary of bounded curvature [31]. 

A Steiner-type formula 

Let A C R d be a non-empty closed set and denote by exo(H) the exoskeleton of A, that is, 
the set of all x £ R d \ A which do not have a unique nearest point in A. Then Theorem 1G 
in [8] says that T~L d (exo(A)) = 0. Thus, % -almost every point x in M. d \ A has a unique 
nearest point in A, denoted by tta(x). The (reduced) normal bundle N(A) C R d x § d_1 
of A is given by 

N(A) := | (ttaOz), :xeR d \(AU exo(H)) J, 

where here and below || • || stands for the usual Euclidean distance and § d_1 stands for 
the Euclidean unit sphere in R d . Lemma 2.3 in [10] implies that N(A) is a countably 
{d— l)-rectifiable subset of R d x S d_1 in the sense of Federer [7], Paragraph 3.2.14. 

Let A be as above. The reach function of A is a strictly positive function on N(A) 
defined as 

8{A, x, n ) := inf{r > 0 : x + rn £ exo(H)} 
for all (x,n) £ N(A). The reach of A is 

reach(yl) := inf{£(H, x,n): (x,n) £ N(A)} 

with the convention that reach(H) = +oo if S(A,x,n ) = +oo for all (x,n) £ N(A). The 
set A is said to be of positive reach if reach(H) £ (0,-t-oo]. In particular, if A is convex, 
then reach(H) = +oo, and vice versa. We also remark that any compact d-manifold with 
C 2 -smooth boundary has positive reach; cf. [ 10 ]. 

If A* denotes the closure of the complement of A, we see that N(dA) := N(A)L) N(A*) 
and we define the extended normal bundle of A as N e (A) := N(A) U TN(A*), where T 
is the reflection map T : x § d_1 —► R d x § d_1 , (a;,n) (a—n). Further, denote the 
reach function of A in this context by <5 + (7l,-,-) G [0,+oo] and define the interior reach 
function S~(A,x,n) := — 8(A*,x, —n) £ [—oo, 0 ] for (x,n) £ x § d_1 . 

From Theorem 5.2 in [10], we know that for each A as above there exist uniquely 
determined signed measures i/q, ..., Vd-i on R d x § d_1 , the so-called extended support 
measures of A, vanishing outside of N e (A), such that the Steiner-type formula 

p d—1 p pS + (A,x,n) 

/ f(x)dx = y^ uj d -j / / r d ~ J ~ 1 f(x + rn) druj (d(x, n)) (2.3) 

jR d \OA j=Q JN e (A) J8-(A,x,n) 

holds for any non-negative measurable bounded function / : —> R with compact sup¬ 
port. Here, for integers j > 0, lo 3 = jnj := 2 tt^ 2 /T(j/2) stands for the surface content of 
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the j -dimensional unit sphere. The signed measures uq, ..., Vd-i encode in some sense 
the singularities of the boundary of A. Although this is not visible in our notation, we 
emphasize that the measures uq,.... v ( i-\ depend on A. 


Admissible sets 

Following [10], we denote by 

d + A := {x £ dA : (x,n) £ N(A) for some n £ § d_1 } 

the positive boundary of A and define Nor(A,:c) := {n £ § d_1 : (x,n) £ IV(A)} for x £ 
d + A. The normal cone at x £ d + A is then nor(A, x) := {an : a > 0,n £ Nor(A,x)} and 
we put 

d ++ A := {x £ d + A : dim(nor(A, x)) = 1}, (2.4) 

where dim(B) denotes the dimension of the affine hull of a set B C R d . Clearly, d ++ A is 
the disjoint union of d 1 A and d 2 A , where 

d k A = (x £ d ++ A : card(Nor(A, x)) = k}, fc£{l,2}. (2.5) 

Let us recall from [14] that a closed set A C R d is called gentle if: 

(i) 'H d ~ 1 (N e (A) fl (B x S d_1 )) < oo for all bounded Borel sets B C R d , 

(ii) for H d_1 -almost all x £ dA there are non-degenerate balls B\ and B 2 containing x 
and satisfying B\ C A and int(l? 2 ) C \ A, where int(l? 2 ) stands for the interior 
of B 2 . 

These assumptions ensure, for example, that 'H d ~ 1 {dA\d + A) = 0; cf. equation (5) in [14]. 
The positive boundary of any closed subset of M. d is ('H d ~ 1 ,d— 1)-rectifiable [10], and thus 
the boundary of every gentle set is 1)-rectifiable, too. In other words, there 

are Lipschitz maps /,; : R d_1 —> R d ,i = 1,2,... such that 'H d ~ l {dA \ Uj>i /^R^” 1 )) = 0; 
see, for example, [7], Paragraph 3.2.14. In particular, at 'H !i_1 -almost every x £ dA there 
is a unique tangent hyperplane denoted by T x := T x (dA). 

Moreover, we recall from [14] that the extended support measures Vj of gentle sets 
have locally finite total variation measures | v :j \ for all j £ {0,...,d— 1}. In particular, 
\vj\(N e (A)) < 00 if A is compact. 

We now define the class A (d) of admissible sets to be the class of compact sets 4cR d 
which are gentle, regular closed and satisfy r H d ~ 1 (d 2 A) = 0. (Recall that a set is reg¬ 
ular closed if it coincides with the closure of its interior.) Here, the assumption that 
r H d ~ 1 (d 2 A) = 0 simplifies the structure of the measure I'd- 1 , to be exploited later. Regu¬ 
larity excludes sets with lower-dimensional ‘tentacles’ attached (e.g., a ball with attached 
line segments). 

The class of gentle and compact sets is rather general and the support measures Uj of 
such sets simplify to well-known objects in special situations. We introduce the following 
classes of sets: 
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• KL d is the class of convex bodies in R d , that is, compact convex sets A C with 
non-empty interior, 

• lZ d is the convex ring, consisting of finite unions of convex bodies in M . d , 

• Ai d denotes the class of compact d-dimensional manifolds in W 1 with twice differ¬ 
entiable boundary, 

• V d is the family of compact sets Acl d with positive reach having non-empty 
interior, 

• UV d denotes the class of all subsets A = A\ U ■ ■ • U A n c R d , n > 1, for sets 
A ll ..., A n G V d and such that f] i£l Ai G V d for any I C {1,... , n}. 

These classes satisfy the inclusions: lC d C V d , IC d C TZ d , V d C UV d , M d C V d and 
lZ d CUV d . If A G K. d , then the extended support measures Vj are related to the gen¬ 
eralized curvature measures of A considered in convex geometry; cf. [26] . A similar com¬ 
ment applies if A G V d is a set with positive reach, for which curvature measures have 
been introduced in [6]. In both cases, it holds that d + A = dA. If A G K, d then A satisfies 
'H d ^ 1 (<9 2 A) = 0. The set classes KL d and V d only contain gentle sets. For the set classes lZ d 
and UV d , curvature measures are defined by additive extension, while for A4 d curvature 
measures are defined via classical differential-geometric methods; see Section 3 in [10] 
for a detailed discussion. Moreover, for sets A £UV d we have that 'H d ~ 1 (dA\d + A) = 0 
(see [10], page 251). Furthermore, if A G lZ d is regular closed, then A is gentle according 
to [14], Proposition 2. Additionally, many UV d -sets (namely those admitting a so-called 
non-osculating representation) are gentle by Proposition 3 in [14]. 


Admissible score functions 

We next consider the collection S of admissible score functions. By this we mean the 
collection of all real-valued Borel measurable functions £(x,X,dA) defined on triples 
(x, A,ch4), where X C R d is locally finite, x G X, A G A (d), and such that £ is trans¬ 
lation and rotation invariant. By the latter two properties, we respectively mean that 
£(x, X, dA) = f(x + z,X + z,dA + z) and that £(x, X, dA) = ^(i9x,i9X,i9(dA)) for all 2 G 
R d , rotations $GSO(d) and input (x,X,dA). If x ^ X, we abbreviate £(x, X U {x}, dA) 
by £(x,X,8A). 

We recall now the concept of a stabilizing functional which was introduced in [22-24] 
after earlier works [12, 15]; see also the surveys [27, 30]. Roughly speaking, a functional 
stabilizes if its value at a given point only depends on a local random neighbourhood and 
is unaffected by changes in point configurations outside of it. Following [31], we need to 
go a step further in the standard framework to account for the dependency of functionals 
£ G S on surfaces. 

To make this precise, denote by B r (x) the Euclidean ball of radius r G (0,oo) and 
centre x G R d and by a homogeneous Poisson point processes on R d of intensity 

r G (0, 00 ). Say that £ G 5 is homogeneously stabilizing if for all r G (0,oo) and all (d — 
l)-dimensional hyperplanes H , there is an almost surely finite random variable R := 
i?(£,'P]l om , H) depending on £, r p^ om and H, the so-called radius of stabilization, such 
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that 

£(o, P T hom n b r (o), h) = £(o, (p T hom n b r (o)) u a, h) (2.6) 

for all locally finite sets A C B R { 0) c , where 0 stands for the origin in R d . Given (2.6), 
the definition of £ extends to Poisson input on all of R d , that is, 

£(0 ,V^ om ,H) = lim ((0,? T hom nB r (0),fl). 

1 ->-00 

Given A € A (d), say that £ is exponentially stabilizing with respect to the pair {V\ ,dA) 
if for all x £ R d there is a random variable R := R(£,x,V\,dA), also called a radius of 
stabilization, taking values in [0,oo) with probability one, such that 

fx{x,Vx n B x -i/d R {x),dA) = £ x {x, {Vx n B x -i, dR (x)) U A, dA) (2.7) 

for all locally finite A C R d \ B x -i/d R {x), and the tail probability satisfies 

limsup-log sup P[R{f,x,V x ,dA) > t] < 0. 
t—yoo t A>o,xeR d 

Surface-order growth for the sums (2.2) involves finiteness of the integrated score £a(^ + 
r\~ 1 / d n,VxidA) over r€l. Thus, it is natural to require the following condition; see 
[31]. Given A £ A (d) and p £ [1, oo), say that ^ satisfies the pth moment condition with 
respect to dA if there is a bounded integrable function G^ ,p := G^ ,p,dA : R —> R + with 
fZo rd ~ 1 (&’ P ( r )) 1/P dr < oo and such that for all r £ R we have 

sup sup supE|£a(£ + r - A _ly,d n,7 :> A VJz,dA)\ p < G^ ,p (|r|). (2.8) 

2£R d U0 (x,n)eN e (A) A>0 

Given A £ A(d), recall for "H d_1 -almost all x £ dA that T x := T x (dA) is the unique 
hyperplane tangent to dA at x. For x £ dA, we put H x := T 0 {dA — x). The score £ is 
said to be well approximated by Vx input on half-spaces if for all A £ A (d), almost all 
x £ dA, and all w £ R d , we have 

lim E\aw,\ 1/d (Vx-x),\ 1 / d (dA-x))~f(w,\ 1 / d (V x -x),H x )\=0. (2.9) 

A—)-oo 


General theorems giving first- and second-order asymptotics 

The results asserted in Section 1 are consequences of general limit theorems giving ex¬ 
pectation and variance asymptotics for the statistics (2.2). We first describe the general 
theory and then, in Section 4, show how to deduce the assertions of Section 1. The general 
limit theorems given here extend Theorems 1.1 and 1.2 in [31] to the class of admissible 
sets and they yield the first- and second-order asymptotics for statistics of other surfaces, 
as discussed in Remark (iii) below. 
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For a score function £ £ H, we put 

KS,dA)~ [ r Eao+sn^^^-^WdsH^idx), (2.10) 

Jd'AJ -oo 

where n is the unique unit normal at 0 with respect to We now state a general 

result giving expectation asymptotics for sums of score functions. Let C(dA) denote the 
set of functions on R d which are continuous at all points x £ dA. 

Theorem 2.1. Let A £ A {d) and k £ C(dA). Suppose that £ £ 5 is homogeneously stabi¬ 
lizing (2.6), satisfies the moment condition (2.8) for some p £ [l,oo), and is well approx¬ 
imated by V\ input on half-spaces as at (2.9). Then for m £ {1,2}, we have the following 
weak law of large numbers: 

lim \-( d - 1 '>/ d Ht(P x ,dA)=p(Z,dA) inL m . ( 2 . 11 ) 

A—>-oo 


Next, we turn to variance asymptotics and define for x,x r £ R d , r £ (0, oo), and all 
(d — l)-dimensional hyperplanes H, 


■V 


hom 


H):= 


E£(s, V* om U {a/}, H)Z(x’, V* oin U {x},H) 
- Ef (&, T* om , H)Ef(x', P* om , H). 


Moreover, define er 2 (£,dA) by 

a 2 (f,dA) 

:=p(e,dA) ( 2 . 12 ) 

p p poo poo 

+ / / c«(0 + rn,p+sn;P^,R d - 1 )K(x) 2 dsdrdp'H d - 1 (da;). 

Jd 1 A J R d_1 J—oo J —oo 

The following general result gives variance asymptotics for sums of score functions. 

Theorem 2.2. Let A £ A(d) and n £ C(dA). We assume that £ £ H is homogeneously 
stabilizing (2.6), exponentially stabilizing (2.7), satisfies the moment condition (2.8) for 
some p£ (2,oo) and is well approximated by V\ input on half-spaces as at (2.9). Then 

lim A - (d - 1)/d Var [H^(P x ,dA)}= a 2 ^,dA). (2.13) 

A—^oo 

Some of the applications presented in Section 1 require the limit theory for the non- 
re-scaled sums J2 x ev a To state the result in this case, call a score function 

£ homogeneous of order 7 £ R if for all a £ (0, 00 ), 


f(ax, aX, a(dA)) =a' 1 f(x,X,dA). 
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When £ is homogeneous of order 7 , it follows that 

Y, Z(x,'P\,dA) = \-^ d Ht(Vx 1 dA). 


xGV\ 


Homogeneity, together with the distributional identity V K / Xj = k(x) x ! d V\ gives 
n(£,dA)= [ E £(0 + sn, V^ om , H*" 1 ) ds [ k(x ) 1 - 71 d U d ~ x (dx) 

J— oo Jd 1 A 


/ oo 

E£(0 + sn, V^ om , K d_1 ) ds • H d ~ x ( dA ) 

-OO 


(2.14) 


and 


a 2 (£,cM) 


: [ E^ 2 (0 + sn,iP 1 hom ,M d - 1 )ds f n{x) 1 -^ /d 'H d - 1 {&x) 

J —00 Jd 1 A 

/* r-OO pOO 

+ / / c ^(0 + rn,p + sn;7 ? | lom ,R d_1 )ds dr dp 

J R d_1 J — oo J — oo 

x f K {x) 2 -^ /d n d -\dx) 

J d 1 A 

/ oo 

E £ 2 (0 + sn, 7 ?h° m , !) ds . H d -1 ( 5j4 ) 

-OO 


(2.15) 


/* /-OO /*oo 

+ / / / c «(0 + m,p + sn;7 ? 1 hom ,R d - 1 )dsdrdp-'H^) y (5A). 

./R d_1 J — oo J — 00 


Consequently, with p,(£,cM) and <r 2 (£,cM) as in (2.14) and (2.16), respectively, we have 
under the conditions of Theorems 2.1 and 2.2 that 


lim \~(d-i-7)/d Y £(x,Vx,dA) = n{^dA) 

A —>00 z —' 

x&V\ 


in L m for m € {1,2}, and 


lim A-^-i-^l/^Var V £( 2 ,7> A ,3A) = ct 2 (£, dA). 

A—>00 z ' 

xeV\ 


(2.16) 


(2.17) 


Remarks. 

(i) Convergence of random measures. The methods presented here also yield expec¬ 
tation and variance asymptotics for integrals of the empirical measures 

Y ^x(x-‘Px-dA)S x 

xeP\ 
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against elements of C{dA) (here, S x stands for the unit-mass Dirac measure at x). The 
details of this extension are straightforward and may be found in, for example, [30], which 
deals with volume-order asymptotics for sums of score functions. 

(ii) Central limit theorems. Say that £ decays exponentially fast with respect to the 
distance to dA if for all pG [l,oo) the function G^’ p defined at (2.8) satisfies 

limsup |u | -1 logG^’ p (|u|) < 0 . (2.18) 

| -zx | —>-oo 


Let <!>(•) denote the distribution function of a standard normal random variable. If £ G H 
decays exponentially fast as in (2.18) and if £ satisfies the moment condition ( 2 . 8 ) with 
p = 3, then by Theorem 1.3 of [31], the statistics (2.2) satisfy a central limit theorem 


sup 

xGM 


P 


Ht(V Xl dA)-EHt(V Xl dA) 
VVar [H£(V x ,dA)] 


<I>(x) 


<r(A) 


with rate function 


r(A) :=c(logA) 3 d+ 1 A (d - 1 )/d (Var[^(PA,c)A)]) 3/2 , 

where c > 0 is a constant not depending on A. In particular, if er 2 (£,cM) is strictly 
positive, then r(A) = c(log A) 3 d+ 1 A -( - d-1 - ) / 2d . This is the case for the examples in Section 1, 
provided that k = 1 and that dA contains a C 2 -smooth subset with positive P d_1 - 
measure. 

(iii) Further applications of general results. Theorems 2.1 and 2.2 have applications to 
statistics of surfaces going beyond those arising in Poisson-Voronoi approximation. For 
instance, these general theorems provide the limit theory for functionals of surfaces of 
germ-grain models including, for example, the limit theory for the number of exposed 
tangent points to Boolean models, as described in Section 3.2 of [19]. Another application 
of the general theory involves the number of maximal points in a sample, which goes as 
follows. A point x GV X is called maximal if the Minkowski sum (R + ) d ® x contains no 
other point of Pa besides x, that is, if ((R + ) d ® x) HP x — {x}. The number M x of maximal 
points of V x has attracted considerable interest in the literature; see [1-3, 5, 11, 31]. 
These works restrict to domains A that are either piecewise linear, convex or smooth. 
We may use Theorems 2.1 and 2.2 to unify and extend these results to domains A which 
are admissible sets, as illustrated by the following statement, whose proof follows from 
modifications of the proof of Theorem 2.5 in [31] and is left to the reader. Let n be a 
density supported on A := {(u,v) G R d_1 xl:a£D,0<»< /(u)}, where D C R d_1 and 
/ : D —> R, and assume that A is an admissible set, that is, A G A (d). We further assume 
that the partial derivatives of / exist a.e. and are bounded away from zero and infinity. 
If P x is a Poisson point process whose intensity measure has density \n with respect to 
Lebesgue measure then there are constants cn G (0, oo) and C 12 G [0, oo) depending only 
on d , k and A such that 


EM a ~chA 1 - 1/ ' 1 and Var[M A ] ~ ciaA 1 ” 1 ^. 
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3. Proofs of Theorems 2.1 and 2.2 

To keep the paper self-contained, we give three preparatory lemmas pertaining to the 
re-scaled scores £ A , A > 0. These are re-formulations of Lemmas 3.1-3.3 in [31], which we 
adopt to our more general set-up. The following lemmas do not require continuity of k 
but instead use that a.e. x £ is a Lebesgue point of k. that is to say 

— [ \k{v)-K( x)\dy 

£ JB e (x) 

tends to zero as e J. 0. Given x £ d 1 A, with d x A defined at (2.5), recall that H x := 
To(dA — x) is the unique tangent hyperplane to dA — x at 0 with unit normal n{x). Let 
0 X denote a point at the origin of H x . 

Lemma 1 . Fix A £ A (d). Assume that £ is homogeneously stabilizing as at (2.6), sat¬ 
isfies the moment condition (2.8) for some p £ (l,oo) and is well approximated by V\ 
input on half-spaces (2.9). Then for all x £ d 1 A, w £ and r£l w have 

lim E£ A (x + rX~ 1/d n(x) + X~ 1/d w,V\,dA) = E£(0 X + rn(x) + w,V^?^.,H x ). (3.1) 

Lemma 2 . Fix AGA(d). Assume that £ is homogeneously stabilizing as at (2.6), sat¬ 
isfies the moment condition (2.8) for some p£ (2,oo), and is well approximated by V\ 
input on half-spaces (2.9). Given x £ d l A, v £ and r £ R, we put for A £ (0,oo), 

X A := £ A (x + r\~ l/d n{x),Vx U {x + rX~ 1/d n(x) + X~ 1/d v},dA), 

Y\ := £ A (x + rA _1 ^ d n(a:) + X~ 1 ^ d v,V\ U {x + rX~ 1 ^ d n(x)},dA), 

X ■■= £(0 X + rn(x),V*°^ U {0^ + rn(x) + v},H x ), 
y ~ £(0 X +rn(x) +v,V^G{O x +rn(x)},H x ). 

Then lim A ^ 00 E[X A F A ] = E[XY]. 

Lemma 3. Fix A £ A (d). Let £ be exponentially stabilizing as at (2.7) and assume the 
moment condition (2.8) holds for some p£ (2,oo). Then there is a constant C £ (0,oo) 
such that for all w,v £ and X £ (0,oo), we have 

|E£ a (iu,P a U {«; + A _1 ^ d u}, <9 A)£ a (uj + X“ l/d v, V\ U {«;}, dA) 

- E£ a (w,P x , 5A)E£ a (w + X~ l ' d v, Vx , dA )| 

< C(E£ a (u;, Vx U {w + X~ 1/d v},dA) p ) 1/p 

x (E£ a (w + X~ 1/d v,Vx U {u>}, 3A) p ) 1/p exp(—C' _1 ||u||). 
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|E£a(w,Pa U {w + A 1 / d u},9^4)^A(w + A 1 / d v,V\U {w},dA) 

-¥4\(w,Vx,dA)E£x(w + X~ 1/d v,Vx,dA)\ < cG f>p (|r|) 1/p exp(-c _1 |H|). 

Proof. The first asserted inequality follows as in either Lemma 4.2 of [20] or Lemma 4.1 
of [4]. The second assertion follows from the first assertion together with the moment 
condition (2.8). □ 


Given these auxiliary lemmas, we may now prove the general results. 


Proof of Theorem 2.1. To show (2.11), it is enough to show the expectation asymp¬ 
totics 

lim A-^-^E V £x{x,Px,dA) = n&d A) (3.2) 

A— yoo *■—' 

xeV x 

and then follow the method of proof of Theorem 1.1 of [31] to deduce L m -convergence 
for m € {1, 2}. 

To show (3.2), we first apply the Mecke identity [26], Theorem 3.2.5, for Poisson point 
processes to obtain 


x -(d-i)/d E £x{x,Vx,dA) = X-( d ~ 1)/d 
xGV\ 


E£a(:c, Pa, dA)Xn(x) da: 


= X 1 ^ f E£x(x,'Px,dA)n(x)dx; 

JR d 

recall that we write t;x(x,Tx,dA) instead of £,x{x,Vx U {x},dA) if x^Vx- We now use 
the Steiner-type formula (2.3) to re-write the last integral as 



r d J 1 E£x(x + rn,V\,dA)K(x + rn)dri'j(d(x,n)), 


where for fixed (x,n) £N e (A), T(x,n) := [5 (A,x,n),5 + (A,x,n)\. Upon the substitution 
r = X~ 1 / d r', we obtain that A _ ( d_1 P ci E J2xev x £a(z,Pa ,dA) equals 


d-1 


V^A-^ / (r'f^Mxix + X-^ d r'n,Vx,dA) 

, =0 d N e (A) J\ 1 / d T(x,n) 

(3.3) 

x k(x + A 1//f Vn) dr'vj(d(x,n)). 

To simplify the notation, write r for r'. By the moment assumption (2.8) with p= 1, 
we conclude that, for each j £ {0 ,... ,d — 1}, the integrand is bounded by the product 
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\r\ d 3 1 G' €a (|r|)||At|| 00 , implying that 


f f r d 3 x E£a(x + A 1 ^ d rn,'P\,dA)n{x + \ 1 ^ d rn) druj(d(x, n)) 

JN b (A) J\ 1 / d T(x,n) 


< 


N e (A) J — oo 


r d 3 1 G^’ 1 (|r|)||K || 00 dr|i/j|(d(a:,7i)) 


/ OO 

r ^-iG«’i(|r|)dr. 

-OO 


The integral r d ~ 3 ~ 1 G^’ 1 (\r\) dr is finite by assumption. Moreover, ||k||oo < oo by 
assumption and \uj\(N e (A)) < oo since A G A (d). Consequently, taking the limit in (3.3) 
as A —> oo, it follows by the dominated convergence theorem that only the term j = d— 1 
remains: 


lim A" (d - 1)/d E Y Z x {x,V x ,dA) 

A—xx> z —' 




= 2 / / lim E£ a (:e + A 1 / d rn,Vx,dA) 

J N e (A) J- oo A_>0 ° 

x k(x + X~ 1 ^ d rn)l(r G X 1 ^ d T{x,n))drvd-i{d{x,n)). 


(3.4) 


Here, we use the identity <ui = 2 and we also use that limA-^oo X 1 ^ d T(x, n) = (— 00 , 00 ), 
which holds by construction of N e (A), where the exoskeleton has been excluded. By 
continuity of k on dA 1 we have limA^oo k(x + A~ 1 /,£ Vn) = n(x). Finally, consider the 
limit 

lim E£ x (x + X- 1/d rn,Vx,dA). 

A—>oo 

To identify it, we use translation invariance and the definition of and write 


^x(x + X-^ d m,Vx ,dA) 

= £a(0 k + A ~ 1/d rn,Vx -x,dA- x) 

= C(Ox + rn, X 1 / d (V x - x), A 1 ' d {dA - x)). 


The measure Vd-i concentrates, according to the discussion around Proposition 4.1 of 
[10], Section 4, on the subset dA ++ of the boundary dA where the normal cone is one 
dimensional; recall (2.4). Moreover, since A G A (d), the measure I'd -1 in fact concentrates 
on the subset d x A C d ++ A (see (2.5)), that is to say, on points of the boundary having 
a unique normal vector or tangent hyperplane as in the case of a smooth surface. 

Since £ is well approximated by input on half-spaces, Lemma 1 implies for all (x, n ) G 
N e (A) with x G d l A 1 that the expectation of the latter expression converges to 

lim E£( 0 x + rn , X 1/d (Vx - x), X 1/d (dA - x)) = E£( 0 X + rn, V%?S, H x ). 

A->oo v ' 
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Thus, we obtain from (3.4), 

lim A- (d - 1)/d E Y fx(x,V x ,8A) 

A—>oo L ' 

xGVx 

(3.5) 

= 2 f f E^(O x + rn 1 V'^ 1 H x )n{x)dru d _ l {d{x 1 n)). 

JN e (A) J-oo 

Now, we simplify the last integral and show that it coincides with as given 

in (3.2). First, recall that there is a unique unit normal vector n(x) at each x £ d 1 A and 
define a measure p, d -i on N{A) by 

= J / 1 ((x,n(x))£-)n d ~\dx). 

1 Jo 1 a 

Since A £ A (d) it follows by Corollary 2.5 and Proposition 4.1 in [10] that 

Pd- i(‘) = i / l((x,n) G-)H Q (x,n)n d ~ 1 {d{x,n)), 

z Jn{a) 

where Ho(x,n) is a certain function depending on the so-called generalized principal 
curvatures of A; see equations (2.13) and (2.24) in [10]. Next, write 


/ f{x,n)v d - 1 (d(x,n))= f(x,n)p, d -i(d(x,n)) 

Jnja) Jn(A) 

f(x,n)p d - 1 (d(x,n)) 


IN (A) 
+ 


f T(N(A*)) 


f(x,n)/j, d -i(d(x,n)). 


J N (A)nT(N (A*)) 

According to the discussion before Theorem 5.2 in [10], given a measurable function / on 
R d x § d_1 , we can split the integral over N e (A) in (3.5) into three parts. The projection 
map tti : N(A) —» R d ,(x,n) i->- x has Jacobian also given by H 0 (x,n ) for ’H d_1 -almost 
all (x,n) £ N(A); see [10], Section 3. Combining these facts with the area formula [7], 
Paragraph 3.2.3, applied to in each of the three resulting integrals, which can be 
combined to a single integral over d 1 A 1 we find that 


(Os + rn , V^,H x )k{x) dri/ d _i (d(z, n)) 


'N e (A) J — oo 


[ [ Ef(O x + rn(x),V^,H x )K(x)H 0 (x,n(x))H 0 (x,n(x)) 1 drU d : (da;) 

Jd 1 AJ- oo 

/ OO 

E^(O x + rn(x),V^,H x )K(x)dr'H d ~ 1 {dx), 

-OO 


d 1 A J —c 


18 


C. Thale and J.E. Yukich 


where we also have used the explicit representation of the measure fid -1 as well as the 
fact that H d ~ 1 (dA ++ \ d 1 A) = 0, which holds because A £ A(d). Since £ is invariant 
under rotations, we may replace H x by R d_1 and 0^ + rn(x) by 0 + rn to obtain (3.2) 
from (3.4), as desired. □ 


Proof of Theorem 2.2. Applying the Mecke formula for Poisson point processes, we 
get 


where 


X-(d-i)/d Yai [H^ Vxjd A)} = X 1 ^ [ E£ x (x,Vx,dA) 2 K(x)dx 

Js. d 

+ \ 1+1 / d [ f IiK(x)n(w)dwdx, 
J R d J 1 * 


(3.6) 


h := ECaO r, V x U M, dA)£ x (w,P a U {x}, dA) - E£ A (x, , cb4)E£ A (w,P x , dA). 


The proof of Theorem 2.1 shows that the first integral in (3.6) converges to 
[ f°° Eao x + rn(x),V^ ) ,R d - 1 )\(x)drH d - 1 (dx) = fi(e,dA). 

Jd 1 AJ -oo 

To complete the proof, we show that the second integral in (3.6) converges to the quadru¬ 
ple integral in (2.12). We re-write the integral with respect to x according to the gen¬ 
eralized Steiner formula (2.3), using the notation already introduced in the proof of 
Theorem 2.1. Furthermore, for all (x,n) £ N e (A), let H(x,n) denote the hyperplane 
orthogonal to n and containing x. Given ( x,n) £ N e (A), we re-write the integral with 
respect to w as the iterated integral over H(x,n) and R. This gives 

yl+l/d f f i lK ( x ) K ( w )d xc ] L w 

JR d Jm d 

= x / / / rd ~ 1 ~ il2 

j =1 J (x,n)€N e (A) J r£T(x,n) Jv£H(x,n) J s£R 

x k(x + rn)n{(x + rn ) + (v + sn)) dsdr drvj(d{x, n)) 


with I 2 equal to 

E£a(x + rn,Vx U {(x + rn) + (v + sn)}, 9 A)£a((x + rn) + (v + sn),Vx U {x + rn}, dA) 
- E£ A (x + rn, V x , <9A)E£a ((x + rn) + (v + sn) , Vx , dA). 

We change variables by putting s = A _ 1 //c V, r = A ~ 1 / d r l and v = \~ 1 / d v'. This transforms 
the differential A 1+1 / d dsdndri/j(d(x,n)) into 

ds' dv' dr'vj(d(x, n)), 


je{i,...,d-i} 
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and I 2 into I 3 given by 

I 3 ■= E£a(x + A ~ 1 /d r'n,V x U{(x + A ~ 1 /d r'n) + (X~ 1 /d v' + A ~ x/d s'n)},dA) 
x &((x + A ~ 1 /d r'n) + ( X~ 1 /d v' + X~ 1 /d s'n), V x U {x + A ~ 1 /d r'n}, dA) 

- Ef x (x + A"VVn,P a, dA)E£ A ((x + A^Vn) + (A^V + A“^Vn),&4). 

To simplify the notation, we shall write s, r and v for s', r' and v', respectively. Then 
^ 1 + 1 /d f f j lK ^x)n(w)dxdw 

JR d jR d 

= f j f [r d ~^h (3.7) 

j—\ JN e (A) J X 1 / d T(x,n) JH(x,n) JM. 

x k(x + A ~ x ^ d rn)K((x + A ~ 1 ^ d rn) + (X~ x ^ d v + A _ 1 / d sn)) ds du dw,-(d(x, n)). 

By the second part of Lemma 3, the factor |J 3 1 in (3.7) is dominated uniformly in A by 
an integrable function of (x,n) £ N e (A), seK, v £ H(x,n) and r£l. More precisely, 

\h\ < cG^’ p (\r\) 1/p exp(-c“ VlMI 2 + s 2 ), 

where the constant c is independent of all arguments. Thus for each j £ {1, 1}, 

we have 

/ / / 

JN e (A) J\ 1 / d T(x,n) JH(x,n) JR 

x k(x + \~ 1 / d rri)K,((x + X~ 1 ^ d rn) + (X~ 1 ^ d v + A _ 1 / d sn)) dsdvdrz^(d(x, n)) 



x exp(—c 1 \/|M| 2 + s 2 ) ds di;dri/j(d(x, n)) 


< cIMILIFjK-^A)) f j f r d 1 J G ? ’ p (|r|) 1 /p exp(-c 1 \/|M | 2 + s 2 ) dsdvdr. 
J -00 Jrj — 1 Jr 

Notice that \i/j\(N e (A)) and the triple integral are finite by the assumption that A £ A (d) 
and the moment condition (2.8), respectively. As in the proof of Theorem 2.1, we have 
linrA^oo X 1 ' d T(x, n) = (—00,00). Taking the limit, as A — > 00, in (3.7) and applying the 
dominated convergence theorem, we see that only the term j = d—1 remains. By Fubini’s 
theorem and Lemma 2, this gives 

lim \ 1+1 / d [ [ 

00 J R d Jm d 


Iik(x)k(w) dx dw 
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= 2 / / / / c 5 (0 x +rn,v + 1 )k(x)' 2 dsdr dwd-i(d(x,n)). 

JN e (A) JR d ~ 1 J —oo J — oo 

We can now use the same arguments as in the proof of Theorem 2.1 to show that the 
integral reduces to the quadruple integral in (2.12). This yields (2.13), as desired. □ 


4. Proofs of Theorems 1.1—1.5 


We shall deduce Theorems 1.1-1.5 from the general Theorems 2.1 and 2.2. In each case, 
it suffices to express the relevant statistic as a sum of score functions and to show that 
the score function satisfies the conditions of the general theorems. We anticipate that the 
expectation formula ( 2 . 10 ) could be evaluated explicitly for some of the score functions 
described below. The proof of the positivity of the constants appearing in the variance 
expressions is postponed to Section 5. 


Proof of Theorem 1.1. We first prove the asserted results for the volume functional 
V\(A) : with A £ A (d). For locally finite X C x £ X , define the score 


£ (1) (a;,A,<9A) 


Vol(u(x) n A c ), 
— Vol(u(x) nA), 
0 , 


if v(x) fl dA / 0 , ieA, 
if v\x) fl dA / 0 , A c , 
if v(x) fl dA = 0 , 


(4.1) 


where v{x) := v(x,X) is the Voronoi cell of x based on the point configuration X. In view 
of the limits appearing in our main results, we also need to define scores on hyperplanes, 
that is, on R d_1 . We thus put 


£ (1) 0r,A, 


pd-1 


rvoiM^nRfT if v{x) ni 

)“<-VolMaOnMl- 1 ), ifw(s)nl 

[ 0 , if v(x) ni 


d—1 


/0,i€ R_ 

^ 0 , xe 


d—1 


\)d—l 


(4.2) 


= 0 , 


where we recall 
ensure that 


^T 1 


x [ 0 ,oo) and !r 1 := 


x (—oo,0]. These definitions 


Vx(A)-Vol(A)= ]T & 1 Xx,P x ,dA) = A " 1 ]T £\x,P x ,dA), 

x&Vx xGV x 


where we use that is homogenous of order d. We wish to deduce the volume asymp¬ 
totics for V\{A) by applying the limits (2.16) and (2.17) with 7 = d and with £ set to . 
It suffices to show that the score j s homogeneously stabilizing (2.6), exponentially 
stabilizing as at (2.7), satisfies the moment condition (2.8) for p = 1 and some p £ (2, 00 ), 
and is well approximated by V\ input on half-spaces as at (2.9). The first three conditions 
have been established several times in the literature; see the proof of Theorem 2.2 of [31]. 

To show that is well approximated by V\ input on half-spaces as at (2.9), it suffices 
to slightly modify the proof of the analogous result in Theorem 2.2 of [31]. For the sake 
of completeness, we provide the details as follows. 
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By definition of A (d), almost all points of dA belong to d l A and it so suffices to show 
(2.9) for a fixed y € d 1 A. Translating y to the origin, letting V\ denote a Poisson point 
process on R d , letting dA denote dA — y , and using rotation invariance of it is 
enough to show for all w £ R d that 

lim E|fW (w, \ 1/d V x ,\ 1/d dA) - f W {w, X 1/d T x , K d_1 )| = 0, 

A—>-oo 

where R d_1 is the unique hyperplane tangent to dA at the origin. Without loss of gen¬ 
erality, we assume, locally around the origin, that dA C ffiil -1 . Fix e > 0 and w £ R d . We 
note that there is a constant L £ (0, oo) such that 

sup(E[£ (1) («;, \ 1/d V x , X 1/d dA) 2 }) 1/2 < L 
\>o 


and 

sup(E[£W(u;, \ l ' d Vx ^ d ~ 1 ) 2 ]) 1/ ' 2 < L. 

A>0 

Let v(w, X^-^Vx) be the union of v(w,X 1 ^ d V\) and all the Voronoi cells adjacent to 
v(w , X 1 /d Vx ) in the Voronoi mosaic of Vx- For all r £ (l,oo), consider the event 

Ei(X,w,r) := {diam(u(uj, A 1 ^^)) < r}, (4.3) 

where diam(-) stands for the diameter of the argument set. Lemma 2.2 of [18] shows there 
is ro :=ro{e,L) such that for r£ [r 0 ,oo) and A large we have P(Ei(X,w,r) c ) < (e/2 L) 2 . 
It follows by the Cauchy-Schwarz inequality that 

lim E|(fW (w, X 1 /d Vx , X 1 /d dA) - ^ ( 1 ) (w, X 1 /d V x ,M d_ 1 ))l (£’ 1 (A, w, r 0 ) c )| < e. 

A—>-oo 

By the triangle inequality and the arbitrariness of e, it is therefore enough to show that 
lim E| (£« (u;, X 1 /d Vx, X 1 /d dA) - ^ (w, A 1 /d P A , R d_ 1 ))l(^i(A, w, r 0 ))| < e. (4.4) 

A—>-oo 

By the way that y was chosen, 0 is a point in d x A and thus has a unique normal vector. 
We first assume w £ R_ 1 ; the arguments with w £ are nearly identical. Moreover, 
we may assume w £ X 1 ' d A for A large. Consider the (possibly degenerate) solid 

A a (w) := A a (w, r 0 ) := (Rt 1 \ X 1 /d A) n B ro (w). (4.5) 

Recalling that dA is ( r H d ~ 1 ,d— 1) rectifiable, it follows that almost all of dA is contained 
in a union of C 1 sub-manifolds of [7], Theorem 3.2.29. Since 0 is a point of d x A, 
it follows that the maximal ‘height’ h x (w,ro ) of the solid A A (w,ro) with respect to the 
hyperplane R d_1 satisfies lim^oo hx(w, ro) =0 for fixed w and ro (see also the linear 
approximation properties of rectifiable sets summarized in Chapter 15 of [17]). Hence, 

Vol(A x(w,r 0 )) = O(hx(w,r 0 ) ■ rfi- 1 ) 
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and so for large A we have Vo1(Aa(u>, ro)) < e. On the event Ei(X,w,ro), the difference 
of the volumes v(w, X^^Vx) D X 1 / d A c and v(w, is at most Vo1(Aa(w>, ro)). 

Thus, for large A we get 

E\(^ 1 \w,X 1 / d Vx 1 X 1 /d dA) - ^(w^X^Vx^-^niE^w.ro)^ 

< Vol(A A (w;,ro)) < e, 

which gives (2.9) as desired. 

We now prove the asserted results for the surface area functional S\(A). As in [31], 
given X locally finite and an admissible set A C R d , define for x £ X n A the area score 
^^(x, X,dA) to be the 'H d_1 -measure of the (d— l)-dimensional faces of v{x) belonging 
to the boundary of U x exnA v ( x ) there are no such faces or if x £ X fl A, then put 
t;( 2 )(x,X,dA) to be zero). Similarly, for x £ X flRl -1 , put ^ 2 \x, A,R d_1 ) to be the 
■H d_1 -measure of the (d — l)-dimensional faces of v(x) belonging to the boundary of 
Uxe^nR d_1 v ( x )> otherwise ^ (x, X, R d_1 ) is zero. We note that is homogenous of 
order d — 1 and that 

Sx(A)= J2 ^ 2 ) (^Vx,dA). 

xeVx 

We wish to deduce the first- and second-order limit behavior of 5 a (A) by applying the 
limits (2.16) and (2.17) with y = d— 1 and with £ set to £( 2 ). 

It is easy to see and well known that the score £^ 2 ) is homogeneously stabilizing (2.6), 
exponentially stabilizing (2.7), and satisfies the moment condition (2.8) for all p > 1; see, 
for example, the proof of Theorem 2.4 of [31]. To see that is well approximated by 
V\ input on half-spaces (2.9), it suffices to follow the proof of Theorem 2.4 of [31]. For 
sake of completeness, we include the details as follows. 

Fix e > 0 and w £ K rf . By the moment bounds on £^ 2 ) and the Cauchy-Schwarz in¬ 
equality, it is enough to show the following counterpart to (4.4), namely to show that 

lim E\(^(w,X 1 /d Vx,X 1 /d dA)-^(w,X 1 /d rx,^ d ~ 1 ))l(E 1 (X,w,r 0 ))\<C£ 1/2 , (4.6) 

A—>-oo 

where E 1 (X,w,ro) is as at (4.3), and where, as above, the origin is a point of d 1 A — y. 
Define 

E 0 (X,w,r 0 ) ■= {X 1/dr P\ n A\(w,r 0 ) = 0}, 

where A A (iu,ro) is as at (4.5). The intensity measure of X 1 ' d T\ is upper bounded by 
Halloo, yielding for large A that 

P[£'o(A,wJ,r 0 ) c ] < 1 - exp(—H kIIoo Vol(A A ('u>,r 0 ))) < ||k||oo£, (4.7) 

where we used that Vol(A A (w, ro)) < e. 

The two score functions X 1 / d V\, X 1 / d dA) and ^ (w, X^! d V\, K d_1 ) coincide on 

the event Ei(X,w, ro) n£'o(A,w,ro). Indeed, on this event it follows that / is a face of a 
boundary cell of X l t d A\ iff / is a face of a boundary cell of the Poisson-Voronoi mosaic 
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of Rl 1 . (If / is a face of the boundary cell v(w, X^^Vx),™ £ A 1 / d A : then / is also a face 
of v(z, X l / d T\) for some z £ A 1 / d A c . If Xf! d V\ D AA(ro,r 0 ) = 0 , then z must belong to 
Ml -1 , showing that / is face of a boundary cell of the Poisson -Voronoi mosaic of R.1 -1 . 
The reverse implication is shown similarly.) 

On the other hand, since 

E[(^ 2 \w,X lld V x ,X lld dA)-^ 2 \w,X lld P x ,R d - 1 )) 2 l(E 1 (X,w,r 0 ))] = O(l), 

and since by (4.7) we have P[i£o(A, w, r o) c ] < ||/c||oo£, it follows by the Cauchy-Schwarz 
inequality that, as A —> oo, 

E| (£ (2) (w, A 1 /d P A , X 1/d dA) - ^ (w, X l / d V x , R d_ 1 ))l(£u (A, w, r 0 )) | 

= E| (£ (2) (w, X 1/d Vx,X 1/d dA) - e (2) {w, X 1 /d V x , R^" 1 )) 

(4.8) 

x l(£ , i(A,w,r 0 ))l(£ , o(A,w,r 0 ) c )| 

<C , (||k|| 00 £) 1/2 . 

Therefore, (4.6) holds and so f 1 ' 2 ' 1 is well approximated by V\ input on half-spaces as at 
(2.9), as desired. □ 

Proof of Theorem 1.2. Let us first recall that the Poisson-Voronoi mosaic is a nor¬ 
mal mosaic; see [26]. This means that with probability one each £-dimensional face in 
skeb(PVA(^4)) arises as the intersection of exactly d — £ + 1 Voronoi cells. 

Now, given X locally finite, x £ X, and an admissible A C R d , define ^ 3 ,l ^(x,X,dA) 
as 

^{x,X,dA):=-^- Y «*(/) 

f£F e (v (x)) 

/C9(PV a (A)) 

and zero otherwise. Then 

H^(A)= Y ^\x,Vx,dA) = X- e ' d Y £i M) (*,PA ,dA), 

xEV\ x€V\ 

xeA 


where we used that £( 3 A) is homogeneous of order £. We wish to deduce the first- and 
second-order limit behaviour of H^\A) by applying the limits (2.16) and (2.17) with 
7 = £ and with £ set to £( 3 ’^. 

The proof that £^ 3 ’^ is homogeneously stabilizing (2.6), exponentially stabilizing (2.7), 
and satisfies the moment condition (2.8) for all p > 1 follows nearly verbatim the proof 
that has these properties. Indeed the radius of stabilization for £( 3 A) coincides with 
that of £( 2 ). 

To see that £( 3 ’^ is well approximated by V\ input on half-spaces as at (2.9), we 
may follow the proof that £( 2 ) is well approximated by V\ input on half-spaces. Notice 
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that on the event E\{\, w, r 0 ) D Eo(\, w, r 0 ), the scores £( 3,e ')(w,\ 1 / d P\,\ 1 / d dA) and 
^ 3 ’ e '>(w,X 1 / d Vx,M d - 1 ) coincide. As in (4.8), we obtain 

E| (C (3 ’ £) (w, A 1 / d dA) - (w, \ 1/d Vx, M d - 1 ))l(S 1 (A, w, r 0 ))\ 

= E| (C (3,£) (w, X 1/d Vx,X l/d dA) - (w, A 1/d Vx, R d_1 )) 

x l(.Ei(A,w,r 0 ))l(.Eo(A,w,ro) c )| < C , (||k||oo£) 1/2 - 

This gives that ^ 3 X) satisfies (2.9) as desired. □ 

Proof of Theorem 1.3. Given X locally finite, x £ X, and A £ A(d), let us define the 
score (x, X , dA) to be the number of ^-dimensional faces of v(x) := v(x,X) belonging 
to <9 (PVa(A)). Define ^ 4,f ^(x, A,K d_1 ) similarly. Then 

fl(A)= ]T ^{x,Vx,dA). 


We shall show that satisfies the hypotheses of Theorems 2.1 and 2.2, and thus deduce 
Theorem 1.3 from (2.16) and (2.17) with £ set to ^A,t) an d ^ se t to zero (notice that ^A/) 
is homogeneous of order 0). For brevity, write £^ 4 ) for ^AV f or fixed l £ {0,. .., d — 1}. 
Now, ^A) is homogeneously and exponentially stabilizing since its radius of stabilization 
coincides with that for the volume score defined in the proof of Theorem 1.1. The 
number N^(x,V a) of ^-dimensional faces of a Poisson-Voronoi cell v(x) has moments 
of all orders and, therefore, the moment condition (2.8) holds because 

l4V + A- 1/d rn, ^{4,571)1 

< nA)( x + A ~ 1/d rn, Vx U {z})\{v(A l / d x + m, A 1/d P A ) n dA + 0) 

for (x,n) £ N e (A). The expectation of the last factor decays uniformly fast in r, giving 
that £ ( - 4 - ) satisfies the moment condition (2.8) for all p > 1 . 

The arguments in the proof of Theorem 1.1 showing that the surface area score is 
well approximated by Vx input on half-spaces extend to show that ^A) is likewise well 
approximated by Vx input on half-spaces. The guiding idea is that with high probability, 
we have that / is a face of a Voronoi cell v(w) belonging to the Poisson-Voronoi approx¬ 
imation of A 1 / d (A — y ) if and only if it belongs to the Poisson-Voronoi approximation of 
KfJ. -1 . Indeed, this happens on the high probability event that the region ‘between’ the 
boundary of the Poisson-Voronoi approximation of A and R d_1 in the neighbourhood of 
w, must be devoid of points; see the proof of Theorem 1.1. Thus, ^AV satisfies all the 
hypotheses of Theorems 2.1 and 2.2 and this completes the proof of Theorem 1.3. □ 

Proof of Theorem 1.4. Given X locally finite, x £ X, an admissible A C R d , and A o C 
dA, put £ (5/) (x,X,Aq) to be the number of ^-dimensional faces of v(x) if v(x) D A 0 ^ 0 
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^(x,X,A 0 ):=Y,^ 5A (.x,X, A o) 

i =o 


and notice that 

Co a (A 0 ) = ^ ^ 5 \x,V x ,Ao). 

x£_V\ 

We shall show that satisfies the hypotheses of Theorems 2.1 and 2.2, and thus 

deduce Theorem 1.4 from (2.16) and (2.17) with £ set to and 7 set to zero (notice 
that is homogeneous of order 0). The score function is homogeneously stabilizing 
as at (2.6), exponentially stabilizing as at (2.7), and satisfies the moment condition (2.8) 
for all p > 1. This is because each with £ £ {0,..., d— 1} has this property. Also, since 
each is well approximated by V\ input on half-spaces for each £ £ {0 ,,d — 1}, 

it follows that enjoys this property as well. Thus satisfies the hypotheses of 
Theorems 2.1 and 2.2, concluding the proof of Theorem 1.4. □ 


Proof of Theorem 1.5. We start with the iterated volume V A> . Conditioned on PV^ 1 " 1 
the first asymptotic equivalence of Theorem 1.1 yields 

E[pf - P A (1) | pvi 13 ] ~ ciA- 1 / d H d - 1 (a(Pvi 1) )). 

Taking expectations and recalling the equivalence ES a (A) ~ c 2 H d ~ 1 (dA ), we obtain 
EP A (2) - V(A) ~ ciA ~ 1/d {c 2 + 


Next, 

EP A (3) - V(A) = EE[P a (3) - P a (2) I PV^ 2) ] + EE[P a (2) - P A (1) | PV^] + EF a (1) - V{A) 
~ ci A - lld c 2 2 H d -\dA) + ciA - x / d OiH d -\dA) + Cl A- 1 / d H d “ 1 (5A) 
= c 1 c 2 , 2 \- 1/d H d ~ l {dA). 


Recursively continuing this way proves the desired claim, namely 
EV a " ) - V{A) ~ c 1 c 2 ,uX^ 1/d n d -\dA). 

The asymptotic equivalences for ESj^, Eand E/ A ^ follow similarly. 


□ 


5. Variance lower bounds 

We complete the proofs of Theorems 1.1-1.4 by proving positivity of the constants ap¬ 
pearing in the variance expressions. The assumption that dA contains a C-smooth 


26 


C. Thale and J.E. Yukich 


subset with positive (d— l)-dimensional Hausdorff measure is essential for our following 
arguments, but we conjecture that this condition can be relaxed. For example, in [28] the 
author establishes upper and lower bounds on Var[VA(A)] for any compact convex set A 
having non-empty interior, without additional smoothness assumptions. However, it is 
unclear (to us) whether the methods of [28] extend to the more general class of admissi¬ 
ble sets A (d) as well as to the other Poisson-Voronoi statistics considered in Theorems 
1.1-1.4. 

In what follows, we use the standard Landau notation. More precisely, for two functions 
f,g : [0,oo) ->lwe write 

• / = o(g) if for all c £ (0,oo) there exists Ao > 0 such that for all A > Ao, |/(A)| < 

c|s(A)|, 

• / = 0(g) if there exists c £ (0,oo) and Ao > 0 such that for all A > Ao, |/(A) | < c|g(A)|, 
and 

• / = fi(g) if there exists c £ (0, oo) and Ao > 0 such that for all A > Ao, |/(A)| > cg( A). 

Positivity of c 3 and c 4 

Positivity of c 3 is shown in Theorem 2.3 of [31] and it remains to consider c\. For this, 
recall that T C dA is C 2 -smooth, with 'H d ~ 1 (T) £ (0,oo). Recalling A C Q, subdivide Q 
into cubes of edge length 1(A) := ([A 1 ^]) -1 . The number L( A) of cubes having non-empty 
intersection with T satisfies L( A) = H(A^ d_1 ^ d ), as otherwise the cubes would partition 
T into o(A < ' d ~ 1 ^ / ' d ) sets, each of 'H d ~ 1 -measure 0(( A _1 / d ) d_1 ), which when A —> oo gives 
'H d ~ 1 (T) = 0, a contradiction. 

We find a sub-collection Q i,..., Qm of the L( A) cubes such that d(Qi , Qj) > 2\fdl(A) 
for all i,j < M , and M = n(A^ _1 ^ d ), where d(Qi,Qj) stands for the distance between 
Qi and Qj. Rotating and translating 1 < i < M , by a distance at most (Vd/2)l(A), if 
necessary, we obtain a collection Q i,..., Qm of disjoint cubes (with faces not necessarily 
parallel to a coordinate plane) such that 

• d(Qi,Qj) > y/dl( A) for all i,j < M, 

• r contains the centre of each Qi 1 here denoted Xi,l<i< M. 

By the assumed differentiability of T, T n Qi is well approximated locally around each Xi 
by the hyperplane T, := T Xi tangent to T at Xj. By the C^-assumption, the approximation 
is uniform over all 1 < i < M. Making a further rotation of Qi, if necessary, we may 
assume that Ti partitions Qi into congruent rectangular solids. Let Ti coincide with the 
hyperplane IR d_1 . Without loss of generality, we assume dA C R d_1 x (—oo,0], that is, 
dA is ‘beneath’ Ti. 

We now exhibit a configuration of Poisson points V\ which has strictly positive prob¬ 
ability and for which <Sa(A) has variability bounded below by YL(A~ ( ' d ~ 1 ^ d 'H d ~ 1 {T)). Let 
e := e(A) := /(A)/28 and sub-divide each Qi, 1 < i < M , into 28 d sub-cubes of edge length 
e. Sub-cubes within Hausdorff distance 4e of dQi are called ‘boundary’ sub-cubes; if a 
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sub-cube is not ^boundary sub-cube then we call it an interior sub-cube. If each bound¬ 
ary sub-cube in Qi contains a point from V\ , then the geometry of the Voronoi cells with 
centres in Qi and distant more than 4e from dQi is not altered by point configurations 
outside Qi (see, e.g., [21]). 

We assume that X, coincides with the origin and we recall that dA C K d_1 x (— 00 ,0] 
so that points near dA may be parametrized by a pair in M d_1 x (— 00 ,0]). By 2(Z d_1 ) 
we mean the set^of all points in R d_1 having integer coordinates of even parity. Consider 
the sub-cubes Qi having the following properties: 

(a) the boundary sub-cubes each contain at least one point from V\, 

(b) V\ D B e / 100 ((ej, ±e)) consists of a singleton for j £ 2(Z d_1 ), \j\ < 10, or 

(b') V\ D B e / 100 ((ej, e/100)) consists of a singleton for j £ 2(Z d_1 ),|j| < 10 and also 
'P\AB e / 100 ((ej, —e/100)) consists of a singleton for j = 0 and j £ 2(Z d_1 )±l, \j\ < 
10, 

(c) V\ puts no other points in Qi. 

(We remark that the choice of the constants 28 and 100 is arbitrary and that we could 
have used any sufficiently large number.) Events (b) and (b'), which each involve 22 
singletons in 22 small balls, happen with the same probability, which is small but bounded 
away from zero uniformly in A, since k = 1. 

Re-labelling if necessary, let I := {1,..., K} be the indices of cubes Qi having properties 
(a)-(c). It is easily checked that the probability a given Qi , 1 < i < M, satisfies property 
(a) is strictly positive, uniformly in A. This is also true for properties (b)-(c), showing 
that 

Ei£ = n(A (d “ 1)/d ). (5.1) 

Abusing notation, let Q := UiLi Qi and put Q c := [0, l] d \ Q. Let F\ be the er-algebra 
determined by the random set /, the positions of points of V\ in all boundary sub-cubes, 
and the positions of points V\ in Q c . Let Ui, 1 < i < M, be the union of the interior sub¬ 
cubes in Qi. If d = 2, we notice that if (b) happens, then the surface dA\ fl Ui contains 
nearly horizontal edges and the total length of these edges is generously bounded above 
by 30e. Indeed, if (b) happens, the 11 cells centered at the points in V\ fl B e / W0 ((ej, — e)), 
j £ {0, ±2, ±4,..., ±10}, contribute to 9 (PVa(A)) a length roughly bounded by the width 
of Ui plus some negligible corrections. On the other hand, if (b') happens then dA\ fl Ui 
contains 10 sharp peaks, with abscissas roughly equal to {±1,±3,.. .,±9}. In fact, it is 
easily checked that dA\ fl Ui contains at least 18 ‘long’, nearly vertical edges of length 
at least 2e, giving a total edge length of at least 36e. A similar situation holds in higher 
dimensions d> 3. 

Conditional on T\, 'H d ~ 1 (dA\r\Q i ) has variability fl(e 2 ^ -1 ^) = fl(A -2+2 / d ), uniformly 
in i £ I, that is, 

Va,r{H d - 1 (dA x nQ i )\T x } = n(\- 2+2 / d ), i £ I. (5.2) 

By the conditional variance formula, 


Var [S X {A)] = Var[E[S A (A)|.F A ]] ±E[Var[S A (A)|J\]] 
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> E[Var[S A (A)|.F A ]] 

= E[Var [H d -\dA x n 2) + H d -\dA x n 2 C )| F x ]\. 


Given F x , the Poisson-Voronoi mosaic of V\ admits variability only inside Q , that is to 
say, given F x , we have 'H d ~ l (dA x fl 2 C ) is constant. Thus, 


Var[S A (A)] > E[Var [U d ~\dA x n 2 )|.Fa]] 


= E 


Var 


Y j n d ~ 1 {dA x nQ l ) 


L iei 




= E^ Var [U d ~\dA x n QOI-^a], 


iei 


(5.3) 


since, given F x , 'H d 1 (dA x fl Qi),i € I, are independent. By (5.1) and (5.2), we have 
Var[5 A (A)] > cA- 2 + 2 /d E[AT] = Sl(\-( d -V/ d ) 


with some finite constant c€ ( 0 ,oo), concluding the proof that C 4 is positive. 


Positivity of cq and eg 

The general idea is to show that configuration (b') generates a surface which has more 
variability (both in terms of complexity and measure) than the surface generated by 
configuration (b). The details go as follows. For t £ {0,l,...,d— 1} and is I, put 
Si t i := (skeb(PV A (j4))) D Ui, noting that 9 (PV a (A)) fl {/» = Sd-i,i (recall the notation 
introduced in the discussion around equation (5.1)). Let S^i{b) be the ^-dimensional 
skeleton arising from configuration (b) and define Se,i{b') similarly. Henceforth, without 
loss of generality we fix i = 1 and write Si for S(a . Observe that Sd-i(b) consists of a 
single (d — l)-dimensional facet / which is nearly a hypercube of dimension d— 1 (and 
nearly a horizontal edge when d = 2 ). Also, Sd- 2 (b) is the union of 2 (d — 1 ) faces, each 
of which is nearly a hypercube of dimension d — 2 . 

On the other hand, Sd-i(b') is the union of (d — l)-dimensional facets Fj, 1 < j < 
2 (d — 1 ), whose union forms the boundary of a solid hyper-pyramid in R d whose base 
is a translate, up to a negligible perturbation, of Sd-i(b). The boundary of the surface 
Uj=i 11 Fj is °f dimension d— 2 and is the union of 2 (d— 1 ) faces, each of which is nearly 

a hypercube of dimension d— 2. In fact, the boundary of the surface ' 1 Fj is a 

translate, also up to a negligible perturbation, of Sd- 2 (b)' 1 we thus denote the boundary 
of \jf d ~ 1] Fj by Sd- 2 {b'). In other words, we have that 


S d -i(b) = S d - 2 (b) U (int /) 







Statistics of the Poisson-Voronoi approximation 


29 


and 


/2(d—i) 

5 d _ 1 (& , )=^- 2 (6 , )U 1J F j \S d - 2 (b) 

V 3 = 1 

Now, Sd- 2 (b ) and Sd- 2 (b') are indistinguishable from the viewpoint of their combina¬ 
torial complexity, as measured by their lower-dimensional skeletons. Moreover, they are 
nearly indistinguishable from a measure theoretic point of view, since the T^-measure of 
their ^-skeletons nearly coincide (modulo negligible corrections). On the other hand, the 
open facet int / differs significantly from Fj \ Sd- 2 {b) in terms of both combi¬ 

natorial complexity and measure. Indeed, the "H^-measure of the ^-skeleton of the latter 
(facets of a pyramid) is strictly larger than the T^-measure of int / (the base of the pyra¬ 
mid). Likewise, for £ £ {0,1,..., d — 2}, the single facet int / has no ^-dimensional faces, 
whereas \ &d -2 (&) has a non-zero number of ^-dimensional faces. These argu¬ 

ments apply to all skeletons Sd-i,i{b),i £ I. By following nearly verbatim the arguments 
showing that c\ is positive, we get that eg and c$ are positive. 


Positivity of c 10 

We have that Coa(A)) is defined in terms of f^(A),£ £ {0,1,..., d — 1}, and it suffices 
to note that configuration (b') leads to a complexity which is strictly larger than the 
complexity arising from configuration (b). We now follow the arguments that C 4 is strictly 
positive. 
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